Efficiency of tunable band-gap structures for single-photon emission 
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The efficiency of recently proposed single-photon emitting sources based on tunable planar band- 
gap structures is examined. The analysis is based on the study of the total and "radiative" decay 
rates, the expectation value of emitted radiation energy and its collimating cone. It is shown that 
the scheme operating in the frequency range near the defect resonance of a defect band-gap structure 
is more efficient than the one operating near the band edge of a perfect band-gap structure. 
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I. INTRODUCTION 

There has been a fast growing interest in light sources 
that produce individual photons on demand, motivated, 
e.g., by potential applications in secure quantum key dis- 
tribution [ljl and scalable quantum computation with lin- 
ear optics |2| • Such sources would also be valuable in in- 
terferometric quantum nondemolition measurements Q 
and as standards for light detection £| . Different applica- 
tions of course impose different requirements on single- 
photon sources. For example, quantum indistinguisha- 
bility of the photons is not necessary for quantum key 
distribution. On the other hand, photons that are in- 
distinguishable and thus can produce multiphoton inter- 
ferences are needed for quantum computation schemes 
with linear optical elements in general [5j- Any realis- 
tic source will be a result of trade-off among different, 
sometimes conflicting requirements. For example, an im- 
provement of the timing of the generation of the photons 
will broaden their spectral profile. 

Search for realizable single-photon sources has been 
widespread, with the schemes ranging from attenuated 
laser pulses and parametric down-conversion to single 
emitters, and numerous strategies have been explored in 
order to enhance the performance of the sources. Faint 
laser sources provide well-collimated narrow-band out- 
puts but suffer from the Poissonian statistics, i.e., the 
possibility of letting more than one photon being emit- 
ted during a single shot. Parametric down-conversion 
produces a pair of photons at a time, allowing one pho- 
ton to herald the existence of the other one. Due to 
the constraint of phase matching, it is possible to know 
the output trajectory, polarization and wavelength of the 
heralded photon. However, the scheme raises a couple of 
problems. Although the use of pulsed pumps restricts the 
pair emission time to a known interval, the actual mo- 
ment of emission cannot be chosen on demand, because 
the conversion is a random process. Moreover, there is 
always the possibility of producing more than one pair 



*Also at Institute of Physics, National Center for Sciences and 
Technology, 1 Mac Dinh Chi Street, District 1, Ho Chi Minh city, 
Vietnam. 



of photons. To overcome the problem of multiple- photon 
emission while maintaining a high likelihood that one pair 
of photons is generated, a multiplexed system consisting 
of multiple down-converters and detectors has been sug- 
gested, such that it is highly probable that one photon 
pair is created somewhere in the array even at low pump 
intensities [|| . The problem of random production might 
be mitigated with the aid of a storage scheme, where once 
a photon pair has been produced, the heralding photon 
activates a high-speed electrooptic switch that reroutes 
the other photon into a storage loop. The stored pho- 
ton keeps circulating in the loop until being switched 
out at a later time chosen by the user 0- The most 
promising strategy to improve the performance of sin- 
gle emitter sources seems to be enclosing the emitter in 
a microcavity. The Purcell effect can then be exploited 
to shorten the lifetime of the excited emitter, thereby 
narrowing the time window during which the photon is 
emitted. Additionally, if the photon is most probably 
emitted into the fundamental mode of the cavity, which 
is often well-directed in space, unwanted omni-directional 
emission can be suppressed, thereby increasing the out- 
coupling efficiency jj]. 

To go a further step towards a deterministic single pho- 
ton source, enclosing of a single emitter within a tun- 
able planar photonic band-gap structure has been sug- 
gested Q . The transition frequency of the excited emit- 
ter is first inside the band gap near its edge, thus being 
prevented from the spontaneous decay. By exposing the 
solid backbone of the photonic crystal to a short optical 
pulse, the band gap is then shifted, e.g., due to the Kerr 
nonlinearity, in such a way that the emitter transition 
frequency now falls in a region of high density of states. 
The emitter is at once forced to decay rapidly, releasing 
a photon. Unfortunately, the scheme cannot work effi- 
ciently, because the key element, namely, the assumed 
high contrast between the densities of states inside and 
outside the band gap (in the vicinity of the band edge) 
is an artifact of the one-dimensional theory |9( . 

As has been shown recently Q, a substantial contrast 
can be realized, if not a perfect photonic crystal is used 
but a photonic crystal containing a defect, with the work- 
ing frequency range being the range around the defect 
resonance instead of that around the band edge. To gain 
first insight, in Ref. Q only the total decay rate of the 
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excited emitter has been studied. Clearly, this quantity 
cannot be regarded, in general, as a comprehensive mea- 
sure of efficiency of the light source, because the decay 
of the excited emitter can be effectively nonradiative or 
the photon is emitted into an unwanted direction. In this 
paper, we give a more detailed assessment of the perfor- 
mance of the schemes, by investigating the expectation 
value of the emitted radiation energy, its angular distri- 
bution, and the "radiative" decay rate. We shall employ 
a quantum mechanical approach that allows one to fully 
take into account absorption and dispersion of the mate- 
rial media. 

The paper is organized as follows. In Scc.[H]an outline 
of the theoretical background is given and the quantities 
to be calculated are introduced. The primarily numerical 
results are presented in Sec. lIIII and a summary and some 
concluding remarks are given in Sec. II VI 

II. BASIC FORMULAS 

A. Poynting vector and energy flow through a 
surface 

Let us begin with introducing the expectation value 
of the slowly varying Poynting vector of the quantized 
electromagnetic field [HJ, E| 

S(r) = (E ( - } (r) x H (+) (r) - H ( ^(r) x E (+) (r)). (1) 

Here, (• • • ) means the quantum mechanical expectation 
value, and E(+)(r) [E (-) (r)] and HW(r) [H (-) ( r )] are 
the positive (negative) frequency parts of the electric and 
magnetic fields, respectively, 

/>OC 

P(r)= duE(r,w), E^(r) = [E (+) (r)] \ (2) 
Jo 

HW(r)= / dwH(r,w), H (_) (r) = [H (+ >(r)] \ (3) 
Jo 

where 

H(r,w) = E(r,a;). (4) 

Further, let 

E(r, w) = E sourcc (r, w) + E frcc (r, w) (5) 

be the decomposition of the electric field in the free-field 
part and the source-field part, with the free-field part 
being in the vacuum state, 

(...E free (r,o;)> = (Ej ree (r,c i ;)...> = 0. (6) 

In particular, in the far-zone limit y/euir/c — > oo, with e 
being an effectively real (positive) permittivity, the rela- 
tion 

V X Source ( r > W ) = * Ve ~ er X E sourco (r, W ) (?) 



is valid (e r , unit vector in the direction of r), and it is not 
difficult to show that in the far-field zone the relationship 

S(r) = 2ce Vi/(r)e r (8) 

holds, where the expectation value 

/(r) = (E<-)(r)E' + )(r)) (9) 

is commonly referred to as the field intensity. Note that 
Eq. JSJ) is in agreement with Eq. (55) in Ref. jl^ . 

The expectation value of the radiation energy trans- 
ported through an area A can be expressed in terms of 
the Poynting vector as 

W= [ At f daS(r,t). (10) 

J J A 

In particular, when the area is the surface of a sphere in 
the far-field zone, 

da = dc/)d6 sin 6r 2 e r , (11) 
then substitution of Eq. iJSJ into Eq. (fTUf) leads to 

//■7T />27T 
dfiW(fi)=/ ddsinB d<f>W(tt), (12) 
Jo Jo 

where 

/>oo 

W(Q)= dt2ce Q y/er 2 I(r,t) (13) 
Jo 

is the expectation value of the radiation energy emitted 
per unit solid angle |13I | . For a planar geometry, a surface 
of the form of a plane is clearly more appropriate. For a 
plane parallel to the xy-plane, 

sin 

da = d.xdwe 2 = -r 2 d6d6e z (14) 

cost/ 

(e z , unit vector in the positive z direction), from Eqs. © 
and IjlOfl we obtain 

W= d6»sin(9 / d0VF(fi), (15) 

Jo Jo 

where, in contrast to Eq. (|12|) . the upper limit of the 
0-integral is now ir/2. 

B. Spontaneous decay 

Consider a single emitter being placed at position ta 
and having transition frequency lua and (real) transition 
dipole moment c!a- Various experimental realizations of 
single emitters have been considered (see, e.g., Refs. 0,0 
and references therein). Here we do not assert very spe- 
cific assumptions on the emitter except that the emission 
is the result of a spontaneous electric dipole transition, 



3 



which implies that the size of the emitter is small com- 
pared to the radiation wavelength. Very importantly, 
the dipole emitter is surrounded by (linear) macroscopic 
bodies, which can be dispersing and absorbing and of 
arbitrary geometry. 

Basing on quantization of the electromagnetic field in 
causal linear media, one can show that the spontaneous 
decay rate of an electric dipole emitter can be determined 
according to the formula 
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d A ImG(r Aj rA,w A )d 
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(16) 



where G(r, r',u>) is the classical Green tensor of the 
medium-assisted Maxwell-field characterizing the sur- 
rounding media. Note that the atomic transition fre- 
quency wa already incorporates the medium-induced 
level shifts. The intensity of the spontaneously emitted 
light registered at some position r_( not necessarily in the 
far-field zone) and time t reads [l3| 



J(r,i) 



dc 4 



|G(r,r A ,WA)d A | e 



2 -rt 



(17) 



and the Poynting vector can be shown to be (App. [XJ 



S(r,i) 
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Im{[G*(r, r A ,w A )d/ 



x [V x G(r,r A ,w A )d A ]}e 



-rt 



(18) 



Note that transit-time effects have been disregarded in 
Eqs. (JT3 and 02). 

It is worth noting that Eqs. ijTo|) - (|T%)) apply to ar- 
bitrary geometry of the dispersing and absorbing bod- 
ies surrounding the dipole emitter. With the help of 
Eqs. i|l(J|) and l|18|) the radiation energy transported 
through an area located at any distance from the ra- 
diation emitting source can be calculated. In particular, 
when the area of interest is in the far zone, then the en- 
ergy passing the area can be related to the field intensity 
through the relationship JSJl and formulas such as i|12|) 
and l(T3|l can be used. 
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FIG. 1: A single (electric) dipole emitter (position za) em- 
bedded in a planar band-gap structure (a) without defect 
(plate thickness — Ao/4, k = 1, 2, . . . , n — 1) and (b) with 
defect (plate thickness dk = Ao/4 if k = 1, 2, . . . , j — l,j + 



, n — 1 and dj = Ao/2). 



For simplicity we shall identify the Green tensor 
G(r,r',w) in Eqs. ljTtj|l - ((l"%|l with the unperturbed Green 
tensor for a planar multilayer, thereby disregarding local- 
field corrections. That is to say, ^(wa) should be real 
and close to unity [ej(w A ) — 1]- The Green tensor for a 
planar multilayer is well known and can be taken from 
Ref. [l^. Substituting its imaginary part in the coinci- 
dence limit r = r' = r A into Eq. 1161) , we obtain 9] 



C. Spontaneous emission in a planar structure 

Let us consider a single (electric) dipole emitter 
sandwiched between two distributed Bragg reflectors as 
sketched in Fig.^ The Bragg reflectors consist of plates 
of thickness Ao /4 of infinite lateral extension and period- 
ically interchanging low and high complex permittivities 
£~l(<-^) and £h(^) in z direction. The zeroth and nth 
"layers" represent the surrounding of the device. It is 
assumed that the jth layer (1 < j < n — 1) contains the 
emitter. If the thickness of the jth layer is also Ao/4, a 
band-gap structure is obtained; if it is Ao/2 the struc- 
ture can be regarded as a planar photonic crystal with 
a defect or, equivalently, as a Fabry-Perot cavity whose 
mirrors are formed by distributed Bragg reflectors. 



il = jHdA||f(fc„) (19) 



[k|| = (k Xl k y ), in-plane component of the wave vector], 
where 
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with 



^ = yJe S ( WA ) -f, Pj = yjki - fcjf , (21) 
+r 1 e-iWiZK-d,) ± f 3 r ? giftdjj^-l ( 22 ) 



(23) 



and ro =k> A rf A / (3^7teo<t) is the well-known decay rate in 
free space. Further, 



+(-) 



arc the total reflection coeffi- 



cients at the upper (lower) stack of layers, where q = s and 
q=p refer to TM and TE polarized waves, respectively. 

Equation (fTUf) together with Eq. i|2U|) represents the 
total decay rate of the emitter, including the contribu- 
tions from the "nonradiative" decay caused by material 
absorption and from the unwanted laterally guided radi- 
ation. To facilitate the determination of the "radiative" 
decay rate r ra( j that is related to the radiation observable 
outside the device, the bulk and scattering parts of the 
Green tensor are, in contrast to Ref. 0, not separated 
from each other, because both contribute to r ra( j. Note 
that the bulk part of the Green tensor can be singled out 
by setting r\ = r q _ = in Eqs. © and 

Let us assume that the surrounding area of the de- 
vice is vacuum or an effectively nonabsorbing medium, 
i.e., £„(wa) [=£q(wa)] can be regarded as being real. 
Then the radiation that can be collected outside the de- 
vice (though not in their entirety because of the finite 
aperture of the collection optics) results from waves with 
fell < k n (=fc ), because only for these waves j3 n (=/?o) is 
real. Thus, the rate r ra( j related to the emission of these 
waves can be obtained according Eq. I|19l) . by restricting 
the upper limit of the integral to k n (= ko), 



r 



rad 



To 



dfe|| f(fc||). 



(24) 



The ratio r ra( j/r may be regarded, at least for weak ab- 
sorption, as one of the measures of efficiency of the device. 
It is, in a sense, an estimation of the weight of the decay 
channel associated with the emission of radiation that 
can leave the device in principle. However, it comprises 
the radiation on the two sides of the device. 

A quantity that clearly distinguishes between the two 
sides and accounts for the radiation that really escapes 
from the device on the wanted side is the expectation 
value of the radiation energy transported through a dis- 
tant area on this side, as defined by Eq. (|10|) . Let us 
choose a plane z — const above the device in Fig.^as the 
area, so that Eq. (|15|l applies. The relevant Green tensor 
now reads 0] 



G(r,r A ,^ A ) = J d 2 k|, e l ( k '' r '' + ^M(k||), 



(25) 



where 
u(k„) 



1 



(2tt)2 2/3„ 
1 



X;e,e+(k,|)£^(-k||, WA ,« A ), (26) 



(T/3«e fc|| + fe||e z ) = e^„(-k||), 



e ?J k ||) = e fcn x e z = -ef„(-k|| 



£gi(-k||,W A ,2A) 



n/3 



D 



<IJ 



(27) 
(28) 

(29) 



x [ e -.(-k||)e-*^+e+.(-k||)rle*^] 



[ r || = (x,v), r A || =0; £ p = l, 6 = -l; e k v unit vector in k|| 
direction; t* ,^ , generalized transmission coefficients] . For 
real k n and within the stationary phase approximation, 
Eq. approaches |Tt| 

G(r,r A ,Lj A ) = -2nik n ~u('^) — (30) 
r \ r ) r 

as k n r — ► oo. 

Now we can calculate the expectation value of the ra- 
diation energy emitted per unit solid angle, W(Q), by 
substituting in Eq. I|17|l for the Green tensor the far-field 
expression given by Eq. I|30|l and inserting the resulting 
expression for 2"(r, t) in Eq. (|13|) . Integrating W(f2) with 
respect to the polar angle <j> over the (27r)-interval, 



W{9) 



A(f)W{VL), 



(31) 



after some calculation we obtain 



W{6) 



8(r/r 



d 2 Az 2k 2 l sin 2 9 



4 



N 2 



a n I 2 
9p+\ 



d M\(\fr\\„n ,2 , l„« ,2 



4 



N 2 , p ~ 



+ \9s 



(32) 



where 



9' g ± 



n/j 



J0jidj-z A ) 



D a 



[l±r q _e 2l ^ ZA ]. 



(33) 



Note that when the transition dipole moment is normal 
to the area, parallel to the area, or completely random, 
then the relation VK(f2) = W(9)/(2ir) is valid, because 
of the symmetry of the system with respect to (f>. Fi- 
nally, from the quantity sin 8W (9), which is the amount 
of energy emitted per unit azimuthal angle 9, the total 
amount of energy passing an area z = const sufficiently 
far above the device is obtained by further integration 
[cf. Eq. JTSJ], 



r/2 

W = I d9 sm8W(9). 
'o 



(34) 
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From Eq. i'S2\) it is seen that, as expected, only the com- 
ponents of the transition dipole moment in the (xy)-pl&ne 
contribute to the far field observed in the normal direc- 
tion (9 = 0). Note that in free space, where 



W{9) 



3hu>A 



4^2sin 2 + ^i(cos 2 
di di 



9 + 1) 



is valid, Eq. simply yields 

W = i^ A , 



(35) 



(36) 



as it should be. 

So far we have considered the case in which the detec- 
tion area is above the device in Fig.^ Clearly, Eqs. 132|) 
and (|34|l also apply to the case in which the detection 
area is below the device, provided that k n is replaced 
with fco and <7™_(_ is replaced with 
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t q 0/j e il3 ^ A r 

Dai 



l±r q e 2i0i ( dj 



(37) 



Needless to say that 9 is now understood as the angle 
between r and the negative z axis. 

It should be mentioned that the angular distribution 
of the radiation energy of an oscillating classical dipole 
is described by formulas that are quite similar to the for- 
mulas given above and differ from them only in some 
factors. In particular, they have been used to investi- 
gate light scattering in planar three-layer structures 
and to study the efficiency of organic microcavity light 
emitting diode structures ^tJ; with strongly absorbing 
metallic mirrors being taken into account. Since sponta- 
neous decay is basically a pure quantum effect that has 
no classical analog, its consistent description should be 
based on quantum theory, without borrowing from else- 
where. 



III. NUMERICAL RESULTS 

The formulas derived in Sec. Ill CI for the "radiative" 
decay rate and the expectation value of the emitted ra- 
diation energy are valid for any planar (n+ 1) multilayer 
structure for which £j(o;a) — 1 and £„(wa) = £o( w a) is 
real. In the numerical calculations we have set £„(oja) = 
eo(u)^=ej(oj^=ex J (uj) = l and assumed that £h(wa) can 
be modeled by a single-resonance permittivity of Drude- 
Lorentz type according to 



ch(w) = 1 + 



lLO"f 



(38) 



where wp corresponds to the coupling constant, and lot 
and 7 are respectively the transverse resonance frequency 
of the medium and the linewidth of the associated ab- 
sorption line. 

The dependence on the transition frequency of the ex- 
pectation value of the radiation energy emitted into the 
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FIG. 2: The expectation value of the radiation energy emit- 
ted into the half space above the multilayer device in Fig. 
is shown as a function of the transition frequency for a band- 
gap structure (a) without and (b) with a defect layer and 
the orientation of the transition dipole moment (at posi- 
tion za = dj/2) parallel (TV i — > W||) to the layers. The up- 
per stack has 5 periods and the lower stack has 5 (solid 
curve), 6 (dashed curve), and 7 (dotted curve) periods [e„(wa) 
= eo(wa) = £j(wa) = £l(<^a) = 1 and sh(loa) according to 
Eq. with u T = 2Qu , 7 = 10~ 7 cj , and ufp = 1.7299 wr]. 
The inset shows the change of W\\ that results from chang- 
ing lop from ujp = 1.7299 lot [eh (wo) - 4 + i 7.5 x 10~ 10 ] (dot- 
ted curve) to uj p = 1.7529 w T [e H (wo) - 4.0804 + i 7.7 x 10~ 10 ] 
(dashed-dotted curve), with the other parameters being the 
same as above for the dotted curve (b). 
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FIG. 3: The same as in Fig. [2] but for a transition dipole 
moment whose orientation is normal (Wi— ► W±) to the layers. 



half space above the device in Fig. ^ is illustrated in 
Figs. [3 and 13 for a single dipole emitter in a planar band- 
gap structure (a) without and (b) with a defect layer. As 
expected, the radiation energy emitted there responds 
more sensitively to the band-gap structure in the case in 
which the transition dipole moment is parallel to the lay- 
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FIG. 4: The expectation value of the radiation energy emit- 
ted into the half space above the multilayer device in Fig. Qis 
shown as a function of the transition frequency for a band-gap 
structure (a) without and (b) with a defect layer for different 
absorption parameters 7 in Eq. I|380 for £h(w) [7 = 10~ 7 wo 
(solid curve), 7 = 10 _3 o;o (dashed curve), and 7 = 10 -2 cjo 
(dotted curve)]. The upper stack has 5 periods and the lower 
stack has 7 periods, and the transition dipole moment is ori- 
ented parallel to the layers. The other data are the same as 
in Fig. E] The magnification in the inset refers to case (a). 



ers (dA = dA||, W = Wi; Fig- El) than in the case in which 
it is perpendicular to the layers (dA = dA_i_, W — W±; 
Fig- EH- It is further seen that the value of Wn can sub- 
stantially increase, particularly with regard to transition 
frequencies lua > Woi with the number of periods in the 
bottom stack of layers of the band-gap structure contain- 
ing the defect layer. The result clearly shows that highly 
unbalanced band-gap structures with the emitter being 
embedded in an appropriately chosen defect element are 
best suitable for realizing a large amount of radiation 
in the wanted space domain, which greatly helps to col- 
lect an emitted photon. Since photon emission in lat- 
eral directions and absorption of the emitted photon by 
the material cannot be completely suppressed, the upper 
limits W\\ = and W\\ = Hua for a balanced scheme 

and an unbalanced scheme, respectively, are not reached. 
Note that unbalanced band-gap structures with defects 
are actually employed in experimental implementations 
of single-photon sources based on micropost microcavi- 
ties (see, e.g., Ref. 5] and references therein). 

The inset in Fig. [21 points up the effect that when 
the transition frequency is adjacent to the resonance fre- 
quency of the defect layer, then the emission can be 
switched from being (nearly) suppressed to being en- 
hanced, by tuning of the permittivity of the band-gap 
material [here: eh(u>)]. For the data in the figure, the 
appropriate frequency interval is 0.997 loq ... 0.998 luq, 
where the change of W\\ is about of 0.2 Huja ... 0.9 Hoja- 
On the contrary, the change of Wn achievable near the 
band edge of the band-gap structure without the defect 
layer is much less, viz 0.2 Hlua- ■ -0.25 ftiOA at wa- 1.23 wq. 
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FIG. 5: The ratios V l{ /T (solid curve), r|| rad /F (dashed 
curve), and F|| rad /r|| (dotted curve) of the decay rates of a 
dipole embedded in a band-gap structure (a) without and (b) 
with a defect layer are shown as functions of the transition 
frequency. The upper and lower stacks in Fig. have 5 layers 
each, and the transition dipole moment is oriented parallel to 
the layers (r, r ra d ^ Tu , F|| rad ). The other data are the same 
as in Fig. [21 

The result is in full agreement with earlier conclusion 
drawn from an analysis of the total decay rate . 

Figure 21 illustrates the effect of material absorption 
on the emitted radiation energy W\\ . As expected, W\\ is 
seen to decrease with increasing 7. For the data in the 
figure, a noticeable effect is only observed for the band- 
gap structure with the defect layer above the working 
frequency ojq. Clearly, the situation would change if one 
of the mirrors were made of metal, which features a much 
stronger absorption. 

From Fig. [S] it can be seen that, as expected, for a 
balanced device the rate Tn Tg ^ defined by Eq. I|24|) (for 
parallel orientation of the transition dipole moment) is a 
quite reasonable measure of the efficiency (compare the 
corresponding curves in Figs. and [2J . Though its fre- 
quency response is qualitatively the same as that of the 
total decay rate Ty, the values of the two rates can sub- 
stantially differ from each other. It is worth noting that 
the frequency response of r|| ra( j/(2r||) and the frequency 
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response of W\\/(Huja) (solid curves in Fig. 0) agree al- 
most completely for the chosen (small) value of the ab- 
sorption parameter 7. To elucidate this agreement, we 
change in the integral in Eq. I|34|) [together with Eq. (|32[l ] 
the integration variable according to fcy = k n sin#. Thus 
we may write 
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{fin — ~ fey)- Under the assumptions that r* 

r 9 



r 



:r y , 2zA = dj (and real /J,), from Eqs. and (JSSJl it 
follows that 
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h qp j.qgiPjdj |2 ' 



and from Eqs. and it follows that 



Re(±e^C^ 



|1 =p r 9 e i /3j d j |2 



(40) 



(41) 



Comparing Eq. (|23J together with Eq. (J2DJ with Eq. IpEty. 
we see, on taking into account Eqs. jlD|| and (03, that 
the relation 



Trad 

~2T~ 
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hbJA 



is valid, provided that the condition 
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1/3 I 



r«l 2 - 



(42) 



(43) 



is satisfied, i.e., negligibly small absorption within the de- 
vice. When the structure becomes unbalanced, then the 
decay rates T and r rac i change very little (not shown), 
while the emitted radiation energy W can change signif- 
icantly (cf. Fig.0). 

In Figs.[5]-|SJ the emitter is placed in the middle of the 
jth layer of the multilayer device in Fig.^ Figure^illus- 
trates the variation of the decay rates with the position 
of the emitter. For the structure without a defect layer 
[Fig.EJa)], the total decay rate T|| is seen to be minimal 
at the center of the layer and to continuously increase 
as the emitter position approaches an edge of the layer. 
The increase of the total decay rate near the edges ob- 
viously results, for small material absorption, from the 
generation of guided waves that do not contribute to the 
radiation emitted into the half space above the device. 
Accordingly, the dependence on the emitter position of 
the rate F|| ra( j is relatively weak, with the maximum be- 
ing at the center of the layer. For the structure with a 
defect layer [Fig. EJb)], on the contrary, the total decay 
rate Tu exhibits a well pronounced absolute maximum at 




FIG. 6: The ratios ry/ro (solid curve), r|| rad /ro (dashed 
curve), and r|| rad /r|| (dotted curve) of the decay rates of a 
dipole embedded in a band-gap structure (a) without and 
(b) with a defect layer are shown as functions of the dipole 
position za for the dipole transition frequency (a) lua — 1.25 ujo 
and (b) tJA = 1.01 loq. The other data are the same as in Fig. El 



the center of the layer and so does the rate F|| rad . Need- 
less to say that this peak can be advantageously used to 
narrow the time window during which the emission takes 
place. In both cases, the middle of the layer is the op- 
timal position of the emitter. In unbalanced structures, 
though the curves (not shown) are not symmetric with 
respect to the center of the layer, the main features dis- 
cussed above are preserved. 

A large portion of radiated emitted into the half space 
above (or below) the device in Fig. [21 is only one condi- 
tion of efficiency. The emission should take place within 
a solid angle as small as possible, so that a photon can 
be easily collected. A measure of the angular distribution 
of the emitted radiation is the c/)-averaged energy W(0) 
defined by Eq. I|32l) . Its dependence on the transition fre- 
quency within the working range is illustrated in Fig. [7| 
for the band-gap structure (a) without and (b) with a 
defect layer [W{6) h-> W {1 (6»)]. Note that a peak of W {1 (6) 
at 9 = corresponds to a single radiation lobe whereas 
a peak at 7^ corresponds to radiation in the shape of 



8 



W\\(0)/(Hua) 



(a) 



1.22 



0.996 




IV. CONCLUDING REMARKS 

In conclusion, we have given an improved analysis of 
the efficiency of two tunable multilayer schemes recently 
proposed for single-photon emission. In particular, we 
have studied (i) the total decay rate and (ii) the "ra- 
diative" decay rate of the excited state of the emitter, 
(iii) the expectation value of the emitted radiation en- 
ergy, and (iv) the collimating cone of the radiation en- 
ergy. The results clearly show that the scheme operating 
near the defect resonance of a band-gap structure with 
a defect is more advantageous than that operating near 
the band edge of a perfect band-gap structure. 

Throughout the calculations we have assumed that the 
permittivity of the layer containing the emitter is effec- 
tively unity. This restriction may be abandoned by tak- 
ing into account appropriate local-field corrections. For 
an atom embedded at the center of a material sphere of 
permittivity e(u>), the real-cavity model can be shown 
Il8l to lead to the corrected Green tensor in the form of 



Gioc(ta, r A ,WA) = G lo u c (r A , r A , wa) 
n 2 

G refl (r A ,r A ,w A ) (44) 



3e(w A ) 



2e{uo A ) + 1 



for equal positions and 



Gioc(r,r A , wa) 



2e{io A ) + 1 



G(r,r A , w A ) 



(45) 



FIG. 7: The angular distribution of the expectation value 
of the radiation energy emitted into the half space above 
the multilayer device in Fig. Q is shown as a function of 
the transition frequency for a band-gap structure (a) with- 
out and (b) with a defect layer. The upper and lower 
stacks in Fig. Q have 5 layers each [uip = 1.7529 lot ~> 
e H (wo) - 4.0804 + i 7.7 x 10~ 10 ]. The other data are the same 
as in Fig. 



a cone of apex angle of 28 in z direction. For the mul- 
tilayer structure without a defect layer, Fig. |7Ja) shows 
that, for the data used, the radiation is most collimated 
near loa — 1.23 wo- At this transition frequency however, 
the radiation energy emitted in the whole half space has 
not reached its maximum yet (cf . Fig. . The maximum 
is reached at w A ~ 1.24 wo, but there the emission has al- 
ready spread. Although the multilayer structure with a 
defect shows qualitatively the same behavior, the emis- 
sion pattern is much better collimated in this case, as 
can be seen from Fig. Efb) . When the structures become 
unbalanced (not shown), the peak heights are changed 
but their positions remain unchanged. 



for different positions. Here, 



G(r, r', u) = G bulk (r, r', w) + G scat (r, r', to) 



(46) 



is the uncorrected Green tensor, with G bulk (r, r', u>) and 
G scat (r, r', oj), respectively, being the bulk and scatter- 
ing parts, and Gj^" 4 is the corrected bulk Green ten- 
sor j^, 0, ligj . The guess has been made that 
Eqs. H44fl and (|45|l might also hold for other geome- 
tries. If this be right, then inclusion in the analy- 
sis of local-field corrections would be straightforward. 
In particular for nonabsorbing material, i.e., effectively 
real e(cj A ), Eq. ((TBI) with Gi oc (r A , r A , w A ) instead of 
G(r A ,r A ,a; A ) implies that both T and r ra( j should are 
corrected by a factor of {3e(u;a)/[2£(w A ) + l]} 2 , result- 
ing in an unaffected ratio r ra( j/r. Similarly, Eqs. (|13fl and 
(|T7Jl with Gi oc (r A ,r A ,cj A ) and Gi oc (r, r A , u; A ) instead of 
G(r A ,r A ,w A ) and G(r,r A ,o; A ), respectively, imply that 
that VF(J7) and W remain unchanged under the local- 
field corrections. 
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APPENDIX A: DERIVATION OF EQ. (T5J| 

To derive Eq. 118(1 , we begin with the Hamiltonian 01 
H = J d 3 r J dwhLJ&(T,w)f(r,w) 



. (0).+ - 
nuj A a' a 



( 7 t d A EW(r A )+H.c. 



(Al) 



Here, the bosonic fields f(r,o>) represent the dynamical 
variables of the system composed of the electromagnetic 
field and a dispersing and absorbing dielectric medium, 
a = \l)(u\ and <jf = \u)(l\ are the Pauli operators of the 
two-level atom, where |^) is the lower state whose energy 
is set equal to zero and \u) is the upper state of energy 
hui^ ', and d A = (/|d A |u) = (u|d A |^) is the transition dipole 
moment. The (positive frequency part of the) medium- 
assisted electric field in terms of the f(r, uS) reads 

/ fe poo 2 r 

E(+)( r ) = iJJL / dcu % / dV vW(r» 

V 7T£0 Jo C ■/ 

x G(r,r',w)f(r (A2) 

The state vector of the total system at time t may be 
written as [|l(r,w)) = f + (r, w)|{0})] 

+ J d 3 r jTliwe-^CKr.w.tJIlCr.w))!?), (A3) 



where wa = — 5uj, with 5lj being the environment- 
induced frequency shift. By formal solution with respect 
to C;(r,u;,£) of the Schrodinger equation for \i/>(t)), one 
can express C;(r, ui,t) in terms of C u {t') at earlier times 
{t 1 < t) . Making the Markov approximation via the re- 
placement C u (t') — > C u (t), after some straightforward 
calculation one derives 



EW(r)M*)> = e- MAt C u (t) 



/•oo a 

x / dcv — ImG(r,r A ,o;)d A C(^A-^)|{0})|0 (A4) 
Jo c 



[C(x) = 7rJ(a;) + iV/x]. Exploiting the analytical proper- 
ties of the Green tensor and performing the w-integral, 
one arrives at, on neglecting a small nonresonant contri- 
bution, 



EW(r)|V(*)) = e - MAt a t (i)^G(r,r A , WA )d A |{0})|0. 

(A5) 

Substitution of Eq. (|A5|l [C„(i) = e~ r */ 2 ] together with 
Eqs. ©-© into Eq. JTJ eventually leads to Eq. JTHJ). 
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